Spike synchronization of a chaotic array as a phase transition 
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We study how a coupled array of spiking chaotic systems synchronizes to an external driving in a short time. 
Synchronization means spike separation at adjacent sites much shorter than the average inter-spike interval; a 
local lack of synchronization is called a defect. The system displays sudden spontaneous defect disappearance 
at a critical coupling strength. Below critical coupling, the system reaches order at a definite amplitude of an 
external input; this order persists for a fixed time slot. Thus, the array behaves as an excitable system, even 
though the single element lacks such a property. The above features provide a dynamical explanation of feature 
binding in perceptual tasks. 

PACS numbers: 



Temporal versus rate coding for the neural-based infor- 
mation has been open to debate in the neuroscience litera- 
ture fT]. The electrical activity of a single neuron is mea- 
sured by micro-electrodes inserted in the cortical tissue of an- 
imals M,^- This activity consists of trains of action potentials 
or "spikes". In rate coding only the mean frequency of spikes 
over a time interval matters, thus requiring a suitable count- 
ing interval, which seems unfit for fast decision tasks. Tem- 
poral coding assigns importance to the precise timing and co- 
ordination of spikes. A special type of temporal coding is syn- 
chrony, whereby information is encoded by the synchronous 
firing of spikes of selected neurons in a cortical module. 

Let an animal be exposed to a visual field containing two 
separate objects. Since each receptive field isolates a specific 
detail, one should expect a corresponding large set of different 
responses. On the contrary, all the cortical neurons whose 
receptive fields are pointing to the same object synchronize 
their spikes, and as a consequence the visual cortex organizes 
into separate neuron groups oscillatingon distinct spike trains 
for different objects (feature binding) j^. Indirect evidence of 
synchronization has been reached for human subjects as well, 
by processing the EEG (electro-encephalo-gram) data [4]. 

Since feature binding results from the readjustment of the 
temporal positions of the spikes, a plausible explanation is 
based on the mutual synchronization properties of chaotic os- 
cillators; in fact, an erratic distribution of the spike occurrence 
prior to an applied stimulus seems mandatory for adapting to 
different rhythms (|5], and references therein). 

In this Letter we show that spiking chaotic dynamics pro- 
vides a suitable model of the feature binding process, where 
fast synchronization is necessary for completion of a percep- 
tual task [2, 4]. Indeed, the conscious perception of a feature 
(e.g. form or motion) is believed to be due to the synchro- 
nization of the neurons of a cortical module within a limited 
time slot [6]; the time limitation is necessary to leave room to 
successive perceptions. Feature binding implies two temporal 
aspects f?, '4, 's'l, namely, i) it is associated with gamma band 
EEG oscillations (25 msec mean interspike separation), and ii) 
it lasts a few hundred milliseconds. Since we intend to model 
the neurons of a cortical module during a perceptual task, we 



study the transient dynamics of defects in the coupled system. 
A defect is defined as the lack of synchronization of two ad- 
jacent sites. Defect disappearance marks the global synchro- 
nization of a modulus. The two main results of this paper are: 
i) above a critical coupling a chaotic array fully synchronizes 
within tens of interspike intervals; ii) in presence of an input, 
the array displays a collective recognition which lasts for a 
fixed time. 

We focus on a model exhibiting homoclinic chaos (HC) 101 
and cornpare it with a generic chaotic system as e.g. 
Rossler ||81]. Let us consider an array of HC systems with near- 
est neighbor bidirectional coupling ruled by 

±\ = kox\{x2 - 1 - fci sin^ x^) 

±2 = -712^2 ~ 2kox\x2 + gxl +x\+p 

^3 = -714 + 52^2 + 4 + P 

±4 = -J2x\ + zx\ + gx\ + zp (1) 

x\ — -722:5 + zx\ + gx\ + zy 

4 = -/3[2;^-6o+r(/(2;l) + eK-i+a;l+i-2r?Xi)))] 

where ](x\) — ^^^^i and 77* (t) is a variable obeying the fil- 
ter equation 77' = ^d{if ^ ^i) with d — 10^"^. The index i 
denotes the zth site position for i — 1, M. The values of 
parameters are; fco = 28.5714, ki = 4.5556, 71 = 10.0643, 
72 = 1.0643, g = 0.05, p = 0.016, 2 = 10, ^ = 0.4286, 
a = 32.8767, r = 160 and bo = 0.1032. Their physi- 
cal meaning has been discussed in Ref. |7]. This parameter 
choice is by no means critical, since several successive finite 
chaotic windows are available. The mutual coupling consists 
of adding to the xq equation on each site a function of the 
intensity xi (action potential) of the neighboring oscillators. 
Chaos due to the homoclinic return to a saddle focus implies 
a high sensitivity to an external perturbation in the neighbor- 
hood of the saddle |5]. 

At each pseudo-period, or inter spike interval = I SI, HC 
yields the alternation of a regular large spike and a small 
chaotic background. The chaotic background is the sensi- 
tive region where the activation from the neighbors occurs, 
while the spike provides a suitable signal to activate the cou- 
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FIG. 1: Space-time plots of (a) HC systems for e — 0.1 where time 
is in {I SI) units and (b) Rossler systems for e = 0.01 where time 
is in T units. The degree of synchronization is characterized by the 
number of defects Nd estimated far beyond the initial transient: (c) 
100 coupled HC; (d) 40 coupled Rossler. 



pling. We use the generic attribution of HC for a large class 
of systems, with a saddle focus instability 1^. This includes 
a class B laser with feedback |01, a modified Hodgkin-Huxley 
model for action potentials in a neuron membrane fl^ and 
the Hindmarsh-Rose model of generation of spike bursts ifTlll . 
In view of this high sensitivity, we expect that they synchro- 
nize not only under an external driving |12], but also for a 
convenient mutual coupling strength. A quantitative indica- 
tor of this sensitivity is represented by the so-called propen- 
sity to synchronization [13]. Furthermore, for an array of 
coupled systems, the onset of a collective synchronization 
was explored numerically for different values of the coupling 
strength e, showing that the lack of synchronization manifests 
itself as phase slips, that is one spike less or more compared to 
the adjacent site over the same time interval. In a space-time 
plot, space denoting the site position and time the point-like 
occurrence of spikes, phase slips appear as dislocations lll4ll . 
Space-wise, a tiny variation of e changes dramatically the size 
of the synchronized domain fl^. Here we explore the dy- 
namical mechanisms underlying the onset of synchronization, 
both the spontaneous one, occuiTing in a coupled aiTay above 
a critical e in the absence of an external input, as well as the 
stimulated one (semantic response) induced by an external 
stimulus localized at one site. Stimulated array synchroniza- 
tion occurs for couplings below Cc- 

We analyze the interaction of neurons by taking the large 
amplitude spikes as 1 and the remaining chaotic and refrac- 
tory background as 0. Then we define the response time tr as 
the time difference between spike occurrence at neighboring 
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FIG. 2: The degree of synchronization of a coupled array esti- 
mated during the initial transient characterized by entropies 5(0) and 
^(SOO) versus e for (a) 100 coupled HC and (b) 40 coupled Rossler. 



sites starting from the first. When neurons are uncoupled, their 
spikes are uncoiTelated and the distribution of the response 
times (including defects) is spread over a broad range. In- 
creasing the coupling strength correlated clusters appear, be- 
ing however characterized by a large variance in the response 
times and boundary appearance of defects. Indeed, a large 
fluctuation in the response time corresponds to a defect which 
interrupts a sequence of synchronized sites, representing the 
boundary with between two clusters. Increasing further the 
coupling strength we observe the transition to overall synchro- 
nization, accompanied by equal response times for all sites. 
The synchronization is not isochronous, as the response time 
distribution consists of two symmetric non zero peaks. This is 
due to a time lag starting from the two end sites, which obey 
open boundary conditions. In the case of unidirectional cou- 
pling we still have the synchronization without isochronism 
ilSn . but this time with a single t,. peak. 

Spiking systems appear to be the operating units of the 
brain cortex 1 1 ] . In order to find out what is the peculiarity 
of spiking systems in the synchronization processes, we com- 
pare the spiking dynamics of HC with another chaotic system, 
namely Rossler. The oscillations in Rossler are characterized 
by a leading single frequency 1 /T and chaos appears in the 
amplitude of oscillations. Here, no spike but rather phase syn- 
chronization occurs. This implies a time code with a poorer 
resolution, as a spike duration in HC is less than one tenth the 
{I SI), whereas phase is resolved over a sizable fraction of the 
period T. The dynamics is ruled by 
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where a = 0.15, b = 0.2 and c = 10. We analyze the Rossler 
systems in a similar manner: when the oscillation crosses the 
zero line upward we take it as a localized 1 (spike-like); other- 
wise it is 0. Also in this case a non-isochronous synchroniza- 
tion can be observed, however for strong couplings the system 
tend to synchronize isochronously. 

We are interested in the way the defects decay in HC and 
Rossler. In Fig. [T] we report the space-time series starting 
from t — as well as from t = 500 (time being in {I SI) 
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FIG. 3: Entropy S for 100 coupled HC calculated throughout the ini- 
tial time slot, starting from different initial conditions, in the absence 
of external driving and for e = 0.136. 
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FIG. 4: HC array, e — 0.11: (a) Entropy vs time for different ampli- 
tudes of an external signal with period T — 0.5(15'/); (b) long time 
entropy vs input amplitude A, for the same T. 



or T units); space is the linear site sequence. For HC with 
e = 0.1 in both panels of Fig. [T^ we observe the same num- 
ber of defects on average, thus there is no decay of defects. 
We obtain the same results also for longer times (not shown 
here). Defects die out at one site and are born again at another 
one in course of time, due to the temporal relations between 
adjacent sites necessary to induce the escape from the saddle 
region (see the detailed discussion for unidirectional coupling 
in Ref. ifisll ). In Rossler the appearance of defects is a global 
phenomenon: the non zero phase difference between neigh- 
boring sites lasts after a long time (left panel in Fig.[T]5). Once 
the array establishes the synchronization line, defects do not 
reappear again as seen in right panel of Fig. 

We characterize the degree of synchronization in HC and 
Rossler in terms of two quantities, namely, the number of de- 
fects Nd in a space-time interval corresponding respectively 
to the 40 sites and 20 {I SI) and the entropy S. As the con- 
trol parameter we take the coupling strength e. In HC, Nd de- 
creases with e and goes to zero at the critical coupling strength 
Cc = 0.13 (Fig.[T];). In Rossler instead, defects decay mono- 
tonically with the coupling strength (Fig. [TJl) but no sharp 
transition is observed. The degree of order of the system is 
obtained by calculation of the entropy S from the distribu- 
tions of the response times tr in the time series during or far 
beyond the initial transient. The entropy S is defined as 

S{t) = -Y,P{tr,t)\np{U,t) (3) 

tr 

where p{tr, t) is the normalized probability distribution of tr 
evaluated at time t. As done for Nd, we take for p(tr, t) the 
average over intervals of 20 {I SI), starting at t. 

In the case of HC (Fig. |2^) S has a sharp transition at 
e — ec- The transition appears exactly at the critical cou- 
pling for which the number of defects goes to zero (Fig. [T];). 
The transition to synchronization during the transient time has 
similar characteristics as the synchronization far beyond the 
transient (see Fig.|2^). Also in this case the sudden decrease of 
entropy at the critical point is associated with the total disap- 
pearance of defects. When the coupling strength is set above 
Cc, defects decay immediately after the coupling is switched 
on. Indeed, below ec the occurrence of defects hinders syn- 



chronization between sites. The reason is that after the appear- 
ance of a defect, the system needs some time to establish syn- 
chronization again. Such appearance and disappearance of de- 
fects (or equivalently, synchronization and de-synchronization 
of spikes) is responsible for the lack of synchronization in the 
system at small coupling strengths. Finally, we checked that 
the value of the critical coupling is universal for any size of 
the array. The reason for that is the locality of the bidirec- 
tional coupling here considered. 

Now we check if the transition occurs also in Rossler. We 
estimate the entropy S versus e during the initial transient and 
far beyond it, as in the case of HC. In Fig. |2}3 the decrease of 
entropy S versus e is smooth, without sudden changes. It is 
because defects decay steadily during the time evolution. In 
this case, once the synchronization is established, defects do 
not reappear as it happens instead in HC. On the other hand, 
as shown in Fig. |2j3, defects appear during the transient time 
up to large values of the coupling strength. 

Thus, there exists a crucial difference between spike and 
phase synchronization, namely, in HC defects appear and dis- 
appear constantly on average up to ec = 0.13. Beyond Ec 
defects never reappear. On the contrary, in Rossler, there 
are no critical changes in the defect dynamics, and for non- 
zero coupling strengths they always decay in the course of 
the time evolution. The defect dynamics in HC appears as a 
new feature; in fact in sudden symmetry changes of space ex- 
tended Ginzburg-Landau type systems transient defects decay 
always, as shown theoretically |16] and verified experimen- 
tally in He^ [17] and in nonlinear optics [18]. The difference 
in defect dynamics in HC and Rossler highlights the existence 
of the phase transition in the synchronization of the HC sys- 
tems. 

The initial transient is interesting from a neurological point 
of view, since perception implies synchronization within the 
short time slot as remarked above. We must then explore a 
transient phase transition, requiring that all sites be synchro- 
nized within a limited time following the application of an 
external signal. We set the coupling strength slightly above 
the critical value ec = 0.13 and measure the entropy S of the 
initial time slots (< = 0) for different sets of initial condi- 
tions. The results are reported in Fig. [3] There are two things 
to be noticed. First, the way in which the transients differ: 
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FIG. 5: The entropy in time for HC array perturbed at only one site 
by a pulse with amplitude A = 0.09 (left panel) and A = 0.104 
(right panel). In both cases duration of pulse is At = {I SI) and the 
coupling strength e — 0.104. 



the spread of entropies has a maximum some time after the 
coupling has been switched on and before the transition to 
synchronization appears. This deviation is due to the different 
competition histories between the coupled sites. Second, the 
existence of different final values of entropy after the transi- 
tion to synchronization. This difference is determined by the 
existence of different ordered states (attractors) 

From Fig. [3] it appears that the ordered states are reached 
after 20 {I SI). Since feature binding has been related to 
the gamma band neuronal oscillations, with (ISI) ^ 25ms 
iSSl. and a perceptual task requires a few hundred millisec- 
onds 101, such a task has to be over within less than 30{ISI). 

In order to model the perceptual task, the transition to syn- 
chronization has to be accomplished in the presence of the 
external driving. Let us consider a periodic signal applied at 
the first site of the aiTay. We introduce the signal by modula- 
tion of parameter bo as: b — &o(l + Asm2TTt/T). Entropy 
depends on the four parameters e, t, A and T. Leaving to 
Ref. ifioll a detailed study, here we comment on the following 
issues. In Fig.|4^, we plot the time behavior of 5* for e = 0.11, 
below spontaneous synchronization, and show that an input 
with period T = 0.5 {I SI) yields a synchronized array for 
A = 0.2, whereas, for lower A, the interval 20 to 80(ISI) is 
still at high entropy. In Fig.gb, we fix t = 20(/S'/), e = 0.11, 
T = 0.5{ISI) and plot the dependence of S on A. A sharp 
discontinuity occurs around Ac ~ 0.07. Here a novel fea- 
ture emerges: while each individual system is chaotic and 
slightly perturbed by the input, the global array behaves as an 
excitable system, with a sharp transition from a steady high 
value of entropy to a low one. At variance with previous mod- 
els which consist of arrays of excitable individuals, usually 
taken as FitzHugh-Nagumo systems |20], here the excitable 
behaviour emerges as a collective property of the coupled ar- 
ray. To test the excitable property of an aiTay, we show in 
FigShow stimuli of amplitude A (above Ac) and duration At 
induce an excitable collective response. The response has a 
form of synchronization lasting a finite time; then the array 
comes back to its previous state. 

In summary, we have demonstrated a phase transition in 
coupled chaotically spiking systems (HC) and characterized 
this transition in terms of the disappearance of defects at the 



critical coupling strength ec- Below and in the absence of a 
stimulus, defects persist in course of time. On the other hand 
in chaotic non-spiking systems as Rossler, where synchro- 
nization regards the phase, and not the spike positions, defects 
decay at long times even for very small coupling strengths. 
The important feature of the spiking chaotic systems as HC 
is that transition to synchronization in response to an external 
driving is sudden, as it occurs in every day experiences re- 
lated to recalling the memory of events or objects (we recall 
a record in a step manner: either we record it or not). We at- 
tribute this sudden coherent response to a global excitability 
property. 
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